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Mot hodology  ami  <  dig  t  act:  <>r  i  za  t  i  on  studies  during  Fiscal  Years 
1977  and  1976  ( FY  77  and  78)  included  gravity  surveys  in  ten 
valleys  in  Arizona  i  live),  Nevada  (two),  New  Mexico  (two),  and 
California  (one) .  The  gravity  data  were  obtained  for  the  pur¬ 
pose  of  est  imat  i  no  the  gross  st  ructure  and  shape  of  the  basins 
and  the  thickness  of  the  valley  fill.  There  was  also  the  possi¬ 
bility  of  detecting  shallow  rock  in  areas  between  boring  loca¬ 
tions.  Generalized  interpretations  from  these  surveys  were  in¬ 
cluded  in  Frt.ec  Western's  (formerly  I-’ugro  National)  Characteri¬ 
zation  Reports  (PN-TR-26n  through  e). 

During  the  FY  77  surveys,  measurements  were  made  to  form  an 
approximate  1-mile  grid  over  the  study  areas,  and  contour  maps 
showing  interpreted  depth  to  bedrock  were  made.  In.  FY  79,  the 
decision  was  made  to  concentrate  on  verifying  and  refining  suit¬ 
able  area  boundaries.  This  decision  resulted  in  a  reduction 
in  the  gravity  program.  Instead  of  obtaining  gravity  data  on  a 
grid,  the  reduced  program  consisted  of  obtaining  gravity  meas¬ 
urements  along  profiles  across  the  valleys  where  Verification 
studies  were  also  performed. 

The  Defense  Mapping  Agency  (DMA),  St.  Louis,  was  requested 
to  provide  gravity  data  from  their  library  to  supplement  the 
gravity  profiles.  For  Big  Smoky,  Hot  Creek,  and  Big  Sand 
Springs  valleys,  a  sufficient  density  of  library  data  was  avail¬ 
able  to  permit,  construction  of  interpreted  contour  maps  instead 
of  just  two-dimensional  cross  sections. 

In  late  summer  of  FY  79,  supplementary  funds  became  available  to 
begin  data  reduction.  At  that  tine,  inner  zone  terrain  correc¬ 
tions  were  begun  on  the  library  data  and  the  profiles  from  Biq 
Smoky  Valley,  Nevada,  and  Butler  and  La  Posa  valleys,  Arizona. 
The  profile  data  from  Whirlwind,  Hamlin,  Snake  Fast,  White 
River,  Garden,  and  Coal  valleys,  Nevada,  became  available  from 
the  field  in  early  October  1979. 

A  continuation  of  gravity  interpretations  has  been  incorporated 
into  the  FY  8u-81  program,  and  the  results  are  being  summarized 
in  a  series  of  valley  reports.  Reports  covering  Nevada-ntah 
gravity  studies  are  numbered  "K-TR-33-"  followed  by  the  abbrevi¬ 
ation  for  the  subject  valley.  In  addition,  more  detailed 
reports  of  the  results  of  FY  77  surveys  in  Dry  Lake  and  Ralston 
valleys.,  Nevada,  were  prepared.  Verification  studies  were 
continued  in  FY  80  and  81,  and  gravity  studies  were  included  in 
the  program.  DMA  continued  to  obtain  the  field  measurements, 
and  there  was  a  return  to  the  grid  pattern.  The  interpretation 
of  the  grid  data  allows,  the  production  of  contour  maps  which 
are  valuable  in.  the  deep  basin  structural  analysis  needed 
for  computer  modi' 1  i  no  in  the  water  resources  proqram.  Tic 
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gravity  interpretations  will  also  be  useful  in  Nuclear  Hardness 
and  Survivability  (NH&S)  evaluations. 

The  basic  decisions  governing  the  gravity  program  are  made  by 
BMC  following  consultation  with  TRW,  Inc.,  Hr tec  Western,  and 
the  DMA.  Conduct  of  the  gravity  studies  is  a  joint  effort 
between  DMA  and  Ertec  Western.  The  field  work,  including 
planning,  logistics,  surveying,  and  meter  operation  is  done  by 
the  Defense  Mapping  Agency  Hydrographic/Topographic  Center 
(DMAHTC),  headquartered  in  Cheyenne,  Wyoming.  DMAHTC  reduces 
the  data  to  Simple  Bouguer  Anomaly  (see  Section  A1.4,  Appendix 
Al.O).  The  Defense  Mappinq  Agency  Aerospace  Center  ( DMAAC ) ,  St. 
Louis,  Missouri,  calculates  outer  zone  terrain  corrections. 

Ertec  Western  provides  DMA  with  schedules  showing  the  valleys 
with  the  highest  priorities.  Ertec  Western  also  recommended 
locations  for  the  profiles  in  the  FY  79  studies  with  the  pro¬ 
vision  that  they  should  follow  existing  roads  or  trails.  Any 
reauired  inner  zone  terrain  corrections  are  calculated  by  Ertec 
Western  prior  to  making  geologic  interpretations. 
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1.0  INTRODUCTION 

1.1  OBJECTIVE 

Gravity  measurements  were  made  in  Pine  Valley  for  the  purpose  of 
estimating  the  overall  shape  of  the  structural  basin,  the  thick¬ 
ness  ot  alluvial  fil  1,  and  the  location  of  concealed  faults. 
The  estimates  will  be  useful  in  modeling  the  dynamic  response  of 
ground  motion  in  trie  basin  and  in  evaluating  ground-water  re¬ 
sources  . 

1.2  LOCATION 

Pine  Valley  is  located  in  the  southwestern  part  of  Utah  (Fig¬ 
ure  11  in  Beaver  and  Millard  counties.  The  town  of  Milford, 
Utah,  is  approximately  36  miles  (58  km)  east  of  the  valley  on 
Highway  21.  Access  throughout  the  valley  is  good  due  to  an 
extensive  network  of  we  1 1 -ma i ntai ned ,  unpaved  roads.  The  valley 
is  principally  rangeland. 

Pine  Valley  is  bounded  on  the  east  by  the  Wah  Wah  Mountains  and 
on  the  west  by  the  Needle  Range  (Figure  2).  The  area  covered 
by  this  report  lies  between  North  latitudes  38°  00'  and  38° 
45'  and  West  longitudes  113°  30'  and  114°  00'.  The  valley 
is  approximately  12  miles  (19  km)  wide  and  40  miles  (64  km) 
long  . 

1.3  SCOPE  OF  WORK 

A  total  of  439  gravity  stations  was  used  in  this  report.  The 
Defense  Mapping  Agency  Aerospace  Center  (DMAAC)  supplied 
151  gravity  stations  from  its  library,  and  288  new  gravity 
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measurement ;;  were  m.iiio  by  t  ho  Defense  Mapping  Agency  Hydrograph¬ 
ic  Topi  >g  r  aph  l  r  Cciit''i  Geodetic  Survey  Squadron  (  DMAHTC/GSS  )  . 

Pine  Valley  and  Waii  Wah  Valley  were  studied  together,  with  the 
results  presented  in  separate  reports.  The  rectangular  region 
containing  both  valleys  is  the  area  between  North  latitudes  38“ 
0  0'  and  3  8  0  4  5 '  and  We  s  t  longitudes  113°  05'  and  114°  00'. 

There  are  778  gravity  stations  in  the  region.  All  were  used  to 
establish  a  common  regional  gravity  trend  for  the  two  valleys. 

Following  residual  separation,  the  geologic  modeling  of  the  two 
valleys  was  done  independently. 
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2 . 0  GRAVITY  DATA JREDUCT ION 

DMAHTC/GSS  obtained  the  basic  observations  for  the  new  stations 
and  reduced  them  to  Simple  Bouguer  Anomalies  (SBA)  as  described 
in  Appendix  A1.0.  Up  to  three  levels  of  terrain  corrections 
were  applied  to  the  new  stations  to  convert  the  SBA  to  the 
Complete  Bouguer  Anomaly  (CBA).  Only  the  first  two  levels  of 
terrain  corrections  described  below  were  applied  to  the  library 
stations . 

First,  the  DMA AC,  St.  Louis,  Missouri,  used  its  library  of 
digitized  terrain  data  and  a  computer  program  to  calculate 
corrections  cut  to  104  miles  (167  km)  from  each  station.  When 
the  program  could  not  calculate  the  terrain  effects  near  a 
station,  a  ring  template  was  used  to  estimate  the  effect  of 
terrain  within  approximately  3000  feet  (914  m)  of  the  station. 
The  third  level  of  terrain  corrections  was  applied  to  those 
stations  where  10  feet  (3  m)  or  more  of  relief  was  observed 
within  130  feet  (40  m).  In  these  cases,  the  elevation  differ¬ 
ences  were  measured  in  the  field  at  a  distance  of  130  feet 
(40  m)  along  six  directions  from  the  stations.  These  data  were 
used  to  calculate  the  effect  of  the  very  near  relief. 
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3  -  0  GEOLOGIC  SUMMARY 

Fine  Valley  is  a  closed  drai naqe  basin  that  lies  within  the 
Basin  and  Ranqe  Fhys i oq raph i c  Province.  The  rocks  that  crop  out 
in  the  adjacent  mountains  ranqe  in  aqe  from  early  Cambrian  to 
Tert iary . 

Rocks  exposed  in  the  Wah  Wah  Mountains  to  the  east  are  primarily 
nor th- to-northeast  dippinq  early  Cambrian  to  middle  Ordovician 
sedimentary  and  metasedimentary  quartzite  and  shale  and  middle 
Cambrian  to  Permian  limestone  and  dolomite.  Tertiary  mafic  to 
felsic  lava  flows  and  ash-flow  tuffs  crop  out  along  the  south¬ 
eastern  and  southern  edges  of  Pine  Valley.  The  same  type  of 
Tertiary  extrusive  rocks  are  predominant  in  the  Needle  Range  on 
the  west.  Tertiary  intrusive  quar tz-monzon i te  crops  out  locally 
on  both  sides  of  the  valley. 

From  late  Precambrian  to  late  Permian  time,  a  westward  thicken¬ 
ing  wedge  of  clastic  and  carbonate  sediments  was  deposited  in 
western  Utah  alonq  a  north-to-northeast  trending  continental 
shelf.  Thrusting  and  foldinq  began  west  of  this  region  in  the 
Jurassic  and  terminated  to  the  east  with  late  Precambrian  and 
early  Paleozoic  rocks  overthrusting  late  Paleozoic  strata  during 
the  Cretaceous  Sevier  Orogeny  (Thorman  and  Ketner,  1979).  Early 
Tertiary  was  a  time  of  widespread  siliceous  volcanism.  Begin¬ 
ning  in  the  Miocene,  extensional  block  faulting  began  in  western 
Utah.  It  was  accompanied  by  volcanism  that  produced  felsic  and 
ma f i  t —  to- fe 1 s ic  lava  flows. 
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The  present  day  Pine  Valley  is  an  eastward-tilted  graben  whose 
ground-water  basin  is  thought  to  be  interconnected  with  that  of 
Wah  Wah  Valley  to  the  east  (Stephens,  1976).  The  valley  is 
bounded  by  inferred,  north-trending,  high-angle  faults  along  the 
western  side  of  the  Wah  Wah  Mountains  and  eastern  side  of  the 
Needle  Range  (Stephens,  1976). 

The  valley  is  underlain  with  uncemented  to  well-cemented,  older 
Quaternary  basin-fill  deposits  with  interbedded  volcanic  flows 
(Stephens,  1976).  Major  surficial  Quaternary  deposits  in¬ 
clude  fine-grained  lacustrine  and  playa  deposits,  alluvial 
fan-gravels,  eolian  sand,  and  stream-channel  deposits. 
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4.0  I NT  ERP RET AT ION 

The  basis  of  interpretation  is  the  Complete  Bouguer  Anomaly 
( CBA ) .  Drawinq  1  shows  the  CBA  gravity  field  contoured  from 
gridded  values  and  the  location  of  the  gravity  stations. 

Mathematical  treatment  of  irregularly  spaced  data  is  ineffi¬ 
cient.  In  order  to  simplify  the  computer  processing,  the  station 
CBA  and  elevation  data  are  reduced  to  sets  of  values  at  uni¬ 
formly  spaced  points  (nodes)  in  a  geographic  array,  or  grid. 
The  values  at  each  node  are  calculated  from  the  station  data 
within  a  circular  area  around  the  node.  A  bell-shaped  weighting 
function  assigns  greater  weiqht  to  the  nearer  data  points.  The 
node  spacing  is  chosen  to  match  the  average  data  spacing.  A 
1.2-mile  (2-km)  grid  spacing  was  used  for  this  analysis. 

4.1  REGIONAL-RESIDUAL  SEPARATION 

A  fundamental  part  of  the  gravity  interpretation  is  the  separa¬ 
tion  of  regional  effects  from  the  local  effects  of  the  valley 
and  its  fill.  The  CBA  contains  long  wavelength  components  from 
deep  and  broad  geologic  structures  extending  far  beyond  the 
valley.  These  long  wavelength  components,  called  the  regional 
gravity,  were  approximated  by  upward  continuation  of  the 
gravity  field.  Upward  continuations  were  made  to  successively 
higher  elevations  until  the  negative  anomaly  over  the  valley  was 
essentially  smoothed  out.  The  final  continuation  was  calculated 
at  an  elevation  of  140,000  feet  (42,672  m).  This  regional 
field  was  subtracted  from  the  CBA  and  the  resulting  residual 
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gravity  anomaly  war.  adjusted  t.y  a  constant  -9.0  milliqals  so 
that  the  zero  residual  would  t  it  approximately  the  existinq  rock 
c'liit  crops . 

4.2  DENSITY  SELECTION 

The  construction  of  a  aenlouic  m  tel  from  rhe  residual  anomaly 
requires  selection  of  denr.it'.  alues  representative  of  the 
alluvial  fill  and  of  the  underlyinq  rock.  Since  only  very 
generalized  density  information  is  available,  the  geologic 
interpretation  of  the  gravity  data  can  be  only  a  coarse  approxi¬ 
mation.  Seven  borinqs  were  drilled  approximately  160  feet  (49 
m)  into  the  alluvium  during  Verification  studies.  The  average 
of  the  densities  measured  at  the  bottom  of  these  borings  was  2.2 
g /cm-*.  To  account  for  compaction  with  depth  (Woollard,  1962; 
and  Grant  and  West,  1965),  2.3  g/cm^  was  used  in  the  modeling 
process . 

Based  on  the  geology  of  the  surrounding  mountain  ranges,  the 
basement  rocks  underlying  Pine  Valley  are  composed  of  Precam- 
brian  quartzites  and  shales  and  Paleozoic  carbonates  and  sili¬ 
ceous  clastic  strata.  Basement  rocks  throughout  the  Great 
Basin  primarily  comprise  Precambrian  and  Paleozoic  siliceous 
clastic  and  carbonate  strata  with  densities  generally  between 
2.6  to  2.9  q/cm3.  The  Paleozoic  carbonate  rocks  in  Nevada 
and  Utah  are  generally  reported  to  be  relatively  high  in  den¬ 
sity,  on  the  order  of  2.8  g/cm^.  This  value  was  selected  to 
represent  the  density  of  the  basement  rock.  The  density  con¬ 
trast  used  for  modeling  was  -0.50  g/cm^. 
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4.3  MODELING 

Modeling  was  done  with  the  aid  of  a  computer  program  which 
iteratively  calculates  a  three-dimensional  solution  of  gravity 
anomaly  data  (Cordell,  1970).  The  gravity  anomaly  is  represent¬ 
ed  by  discrete  values  on  a  two-dimensional  grid.  The  source 
of  the  anomaly  (the  volume  of  low-density  valley  fill)  is 
represented  by  a  set  of  vertical  prism  elements.  The  tops  of 
the  prisms  lie  in  a  common  horizontal  plane.  The  bottoms  of  the 
prisms  collectively  represent  the  bottom  of  the  valley  fill. 
Each  prism  has  a  cross-sectional  area  equal  to  one  grid  square 
and  a  uniform  density.  A  grid  square  of  1.2  miles  by  1.2  miles 
(2  km  by  2  km)  was  selected  as  representative  of  the  gravity 
station  distribution.  Computations  were  made  for  five  itera¬ 
tions  of  mutually  interactive  prism  adjustments.  The  root-mean- 
square  error  for  the  entire  grid  was  less  than  0.7  milligal. 

The  calculated  thickness  of  the  valley  fill  depends  upon  the 
density  contrast  (i.e.,  fill  density  minus  rock  density)  used. 
Since  neither  density  is  perfectly  known,  nor  even  uniform,  the 
calculated  thickness  should  be  expected  to  contain  a  correspond¬ 
ing  degree  of  uncertainty.  A  source  of  error  in  modeling  Pine 
Valley  as  a  simple  alluvium  basement  rock  system  is  the  wide¬ 
spread  volcanic  material  throughout  the  valley. 

Eight  seismic  refraction  lines  (Table  1)  and  six  borings  (Table 
2)  were  used  as  constraints  in  the  modeling  process.  Their 
locations  are  marked  in  Drawing  2.  The  seismic  refraction 
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lines  located  near  the  mountain  flanks  recorded  high  velocities 
which  may  represent  the  basement  material.  The  alluvial  fill 
material  in  the  center  of  the  valley  is  at  least  1157  to  2006 
feet  (353  to  611  m)  thick  according  to  three  of  the  borings 
described  in  the  literature.  The  thickness  of  basin  fill  (or 
depth  to  rock)  based  on  the  interpretation  of  gravity  data  is 
contoured  in  Drawing  2. 

4.4  DISCUSSION  OF  RESULTS 

The  geologic  structure  of  Pine  Valley  is  interpreted  on  the 
depth- to-rock  contour  map  (Drawing  2).  The  interpretation  is 
based  on  geologic  information  from  published  reports,  analysis 
of  aerial  photographs,  and  geologic  field  reconnaissance  as  well 
as  gravity  data.  The  analysis  of  the  gravity  data  included 
calculation  of  the  second  vertical  derivative  (SVD)  of  the  CBA 
field.  One  property  of  the  SVD  is  that  its  zero  value  marks  the 
steepest  gradients  of  the  input  CBA  field.  This  property  was 
used  to  guide  the  placement  of  faults  in  the  structural  inter¬ 
pretation.  The  interpreted  faults  represent  only  the  major 
fault  systems  which  probably  comprise  many  smaller  fault  zones. 
There  may  be  other  discrete  faults  that  had  a  minor  role  in 
basin  formation,  but  with  displacements  so  small  that  they  were 
not  resolved  by  the  widely  spaced  gravity  data  available  for 
this  study. 

The  depth-to-rock  contours  define  a  major  elongate  north-south 
trending  trough  coincident  with  the  valley  physiography.  These 
contours  appear  to  define  two  subsurface  basins;  a  north¬ 
trending  northern  basin  about  8000  feet  (2438  m)  deep  and  a 
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southern  basin  which  is  about  3500  feet  (1067  m)  deep  and  trends 
sliqhtly  west  of  north  (Drawinq  2). 

The  northern  basin  contains  a  major  qraben.  The  eastern  side  of 
the  graben  has  a  steep  linear  gradient  separating  it  from  the 
Wah  Wah  Mountains.  This  pattern  suggests  there  is  a  continuous 
major  fault  system  along  the  base  of  the  Wah  Wah  Mountains. 
West  of  the  northern  end  of  the  major  graben,  a  complex  pattern 
of  minor  steep  gradients  suggests  an  intricate  arrangement  of 
faults  of  varying  displacements  (Drawing  2,  Figure  3).  Farther 
south,  in  the  vicinity  of  profile  B-B'  (Drawing  2,  Figure  4), 
the  smoother  bedrock  contours  on  the  western  side  of  the  basin 
indicate  that  the  bedrock  dips  gently  eastward  from  the  Needles 
Range  for  several  miles  before  it  is  faulted  down  into  the  major 
graben . 

The  southern  basin  is  more  simple  structurally.  It  is  inter¬ 
preted  to  be  an  eastward  tilted  block  instead  of  a  graben 
(Figure  5). 

The  two  basins  are  separated  by  a  northwesterly  trending  trans¬ 
verse  fault  which  has  allowed  the  northern  basin  to  drop  more 
deeply  than  the  southern  basin.  This  interpretation  is  in 
accord  with  studies  of  surficial  geology  and  geomorphology  which 
indicate  that  displacement  on  the  eastern  basin-bounding  fault 
system  diminishes  southward.  The  fault  appears  to  be  terminated 
before  it  reaches  the  southern  end  of  Pine  Valley  where  no 
evidence  of  major  fault  displacements  is  found  in  Tertiary  lava 
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5 . 0  CONCLUSIONS 

Pine  Valley  gravity  data  indicate  that  the  valley  occupies  the 
down-tilted  portion  of  an  asymmetric  fault  block.  This  fault 
block  is  separated  from  the  Wah  Wah  Mountains  by  a  major  basin- 
bounding  fault  system.  A  narrow  graben,  which  is  about  8000 
feet  (2438  m)  deep,  underlies  the  northern  part  of  the  valley. 

The  calculated  bedrock  depths  are  only  approximations  because 
little  is  known  about  the  actual  density  distribution  in  and 
around  the  valley,  and  the  residual  gravity  anomaly  is  necessar¬ 
ily  based  on  an  interpreted  regional  field.  An  average  density 
contrast  of  -0.50  g/cm3  between  the  alluvium  and  bedrock  was 
used  to  calculate  the  thickness  of  the  valley-fill  material. 
Future  studies  that  acquire  better  density  data  or  measure  actu¬ 
al  depths  to  bedrock  in  deep  parts  of  the  valley  can  be  used  to 
refine  the  gravity  interpretation. 
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A 1  . 0  GENERAL  PR I NC I  PALS  _0F  THE  GRAVITY 
EXPLORATION  METHOD 

A  1.1  GENERAL 

A  gravity  survey  involves  measurement  of  differences  in  the 
gravitational  field  between  various  points  on  the  earth's 
surface.  The  gravitational  field  values  being  measured  are  the 
same  as  those  influencing  all  objects  on  the  surface  of  the 
earth.  They  are  generally  associated  with  the  force  which 
causes  a  1-gm  mass  to  be  accelerated  at  980  cm/sec^.  This 
force  is  normally  referred  to  as  a  1-g  force. 

Even  though  in  many  applications  the  gravitational  field  at 
the  earth's  surface  is  assumed  to  be  constant,  small  but  distin¬ 
guishable  differences  in  gravity  occur  from  point  to  point.  In 
a  gravity  survey,  the  variations  are  measured  in  terms  of  milli- 
gals.  A  milligal  is  equal  to  0.001  cm/sec^  or  0.00000102  g. 
The  differences  in  gravity  are  caused  by  geometrical  effects, 
such  as  differences  in  elevation  and  latitude,  and  by  lateral 
variations  in  density  within  the  earth.  The  lateral  density 
variations  are  a  result  of  changes  in  geologic  conditions.  For 
measurements  at  the  surface  of  the  earth,  the  largest  factor 
influencing  the  pull  of  gravity  is  the  density  of  all  materials 
between  the  center  of  the  earth  and  the  point  of  measurement. 

To  detect  changes  produced  by  differing  geological  conditions, 
it  is  necessary  to  detect  differences  in  the  g r av i ta t ional  field 
as  small  as  a  few  milligals.  To  recognize  changes  due  to 
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geological  conditions,  the  measurements  are  "corrected"  to  ac¬ 
count  for  changes  due  to  differences  in  elevation  and  latitude. 

Given  this  background,  the  basic  concept  of  the  gravitational 
exploration  method,  the  anomaly,  can  be  introduced.  If,  instead 
of  being  an  oblate  spheroid  characterized  by  complex  density 
variations,  the  earth  were  made  up  of  concentric,  homogeneous 
shells,  the  gravitational  field  would  be  the  same  at  all  points 
on  the  surface  of  the  earth.  The  complexities  in  the  earth's 
shape  and  material  distribution  are  the  reason  that  the  pull  of 
gravity  is  not  the  same  from  place  to  place.  A  difference  in 
gravity  between  two  points  which  is  not  caused  by  the  effects  of 
known  geometrical  differences,  such  as  in  elevation,  latitude, 
and  surrounding  terrain,  is  referred  to  as  an  "anomaly." 

An  anomaly  reflects  lateral  differences  in  material  densities. 
The  gravitational  attraction  is  smaller  at  a  place  underlain  by 
relatively  low  density  material  than  it  is  at  a  place  underlain 
by  a  relatively  high  density  material.  The  term  "negative 
gravity  anomaly"  describes  a  situation  in  which  the  pull  of 
gravity  within  a  prescribed  area  is  small  compared  to  the  area 
surrounding  it.  Low-density  alluvial  deposits  in  basins  such  as 
those  in  the  Nevada-Utah  region  produce  negative  gravity  anoma¬ 
lies  in  relation  to  the  gravity  values  in  the  surrounding 
mountains  which  are  formed  by  more  dense  rocks. 

The  objective  of  gravity  exploration  is  to  deduce  the  variations 
in  geologic  conditions  that  produce  the  gravity  anomalies 
identified  during  a  gravity  survey. 
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A  1.2  INSTRUMENTS 

The  sensing  element  of  a  LaCoste  and  Romberg  gravimeter  is  a 
mass  suspended  by  a  zero- lenqth  spring.  Deflections  of  the 
mass  from  a  null  position  are  proportional  to  changes  in  gravi¬ 
tational  attraction.  These  instruments  are  sealed  and  compen¬ 
sated  for  atmospheric  pressure  changes.  They  are  maintained  at 
a  constant  temperature  by  an  internal  heater  element  and  thermo¬ 
stat.  The  absolute  value  of  gravity  is  not  measured  directly  by 
a  gravimeter.  It  measures  relative  values  of  gravity  between 
one  point  and  the  next.  Gravitational  differences  as  small  as 
0.01  milligal  can  be  measured. 


A 1 . 3  FIELD  PROCEDURES 

The  gravimeter  readings  were  calibrated  in  terms  of  absolute 
gravity  by  taking  readings  twice  daily  at  nearby  USGS  gravity 
base  stations.  Gravimeter  readings  fluctuate  because  of  small 
time-related  deviations  due  to  the  effect  of  earth  tides  and 
instrument  drift.  Field  readings  were  corrected  to  account  for 
these  deviations.  The  magnitude  of  the  tidal  correction  was 
calculated  using  an  equation  suggested  by  Goguel  (1954): 

C  =  P  +  Ncos  0  (cos  *5  +  sin  0  )  +  Scos  0  (cos  0  -  sin  0) 
where  C  is  the  tidal  correction  factor,  P,  N,  and  S  are  time- 
related  variables,  and  0  is  the  latitude  of  the  observation 


point.  Tables  giving  the  values  of  P,  N,  and  S  are  published 
annually  by  the  European  Association  of  Exploration  Geophysi- 
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The  meter  drift  correction  was  based  on  readings  taken  at  a 
designated  base  station  at  the  start  and  end  of  each  day.  Any 
difference  between  these  two  readings  after  they  were  corrected 
for  tidal  effects  was  considered  to  have  been  the  result  of 
instrumental  drift.  It  was  assumed  that  this  drift  occurred  at 
a  uniform  rate  between  the  two  readings.  Corrections  for  drift 
were  typically  only  a  few  hundredths  of  a  milligal.  Readings 
corrected  for  tidal  effects  and  instrumental  drift  represented 
the  observed  gravity  at  each  station.  The  observed  gravity 
values  represent  the  total  gravitational  pull  of  the  entire 
earth  at  the  measurement  stations. 

A 1 . 4  DATA  REDUCTION 

Several  corrections  or  reductions  are  made  to  the  observed 
gravity  to  isolate  the  portion  of  the  gravitational  pull  which 
is  due  to  the  crustal  and  near-surface  materials.  The  gravity 
remaining  after  these  reductions  is  called  the  "Bouguer  Anoma¬ 
ly."  Bouguer  Anomaly  values  are  the  basis  for  geologic  inter¬ 
pretation.  To  obtain  the  Bouguer  Anomaly,  the  observed  gravity 
is  adjusted  to  the  value  it  would  have  had  if  it  had  been 
measured  at  the  geoid,  a  theoretically  defined  surface  which 
approximates  the  surface  of  mean  sea  level.  The  difference 
between  the  "adjusted"  observed  gravity  and  the  gravity  at  the 
geoid  calculated  for  a  theoretically  homogeneous  earth  is 
the  Bouguer  Anomaly. 
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Four  separate  reductions,  to  account  for  four  geometrical 
effects,  are  made  to  the  observed  qravity  at  each  station  to 
arrive  at  its  Rouquer  Anomaly  value. 

a.  Free-Air  Effect:  Gravitational  attraction  varies  inversely 
as  the  square  of  the  distance  from  the  center  of  the  earth. 
Thus,  corrections  must  be  applied  for  elevation.  Observed  grav¬ 
ity  levels  are  corrected  for  elevation  using  the  normal  vertical 
gradient  of: 

FA  =  -0.09406  mg/ft  (-0.3086  mil 1 igal s/me ter ) 
where  FA  is  the  free-air  effect  (the  rate  of  change  of  gravity 
with  distance  from  the  center  of  the  earth).  The  free-air 
correction  is  positive  in  sign  since  the  correction  is  opposite 
the  effect. 

b.  Bouguer  Effect:  Like  the  free-air  effect,  the  Bouguer 

effect  is  a  function  of  the  elevation  of  the  station,  but  it 
considers  the  influence  of  a  slab  of  earth  materials  between 
the  observation  point  on  the  surface  of  the  earth  and  the 
corresponding  point  on  the  geoid  (sea  level).  Normal  practice, 
which  is  to  assume  that  the  density  of  the  slab  is  2.67  grams 
per  cubic  centimeter  was  followed  in  these  studies.  The  Bouguer 
correction  (Bc),  which  is  opposite  in  sign  to  the  free-air 
correction,  was  defined  according  to  the  following  formula. 

Bc  =  0.01276  (2.67)  hf  (milligals  per  foot) 

Bc  =  0.04185  (2.67)  hm  (milligals  per  meter) 
where  hf  is  the  height  above  sea  level  in  feet  and  hm  is  the 
height  in  meters. 
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c.  Latitude  Effect:  Points  at  different  latitudes  will  have 
different  "gravities"  for  two  reasons.  The  earth  (and  the 
geoid)  is  spheroidal,  or  flattened  at  the  poles.  Since  points 
at  higher  latitudes  are  closer  to  the  center  of  the  earth  than 
points  near  the  equator,  the  gravity  at  the  higher  latitudes  is 
larqer.  As  the  earth  spins,  the  centrifugal  acceleration 
causes  a  slight  decrease  in  gravity.  At  the  higher  latitudes 
where  the  earth's  radii  are  smaller,  the  centrifugal  accel¬ 
eration  diminishes.  The  gravity  formula  for  the  Geodetic 
Reference  System,  1967,  gives  the  theoretical  value  of  gravity 
at  the  geoid  as  a  function  of  latitude.  It  is: 

g  =  978.0381  (1  +  0.0053204  sin2  0  -  0.0000058  sin220)  gals 
where  g  is  the  theoretical  acceleration  of  gravity  and  is 
the  latitude  in  degrees.  The  positive  term  accounts  for  the 
spheroidal  shape  of  the  earth.  The  negative  term  adjusts  for 
the  centrifugal  acceleration. 

The  previous  two  corrections  (free  air  and  Bouguer)  have  ad¬ 
justed  the  observed  gravity  to  the  value  it  would  have  had  at 
the  geoid  (sea  level).  The  theoretical  value  at  the  geoid  for 
the  latitude  of  the  station  is  then  subtracted  from  the  adjusted 
observed  gravity.  The  remainder  is  called  the  Simple  Bouguer 
Anomaly  (SBA).  Most  of  this  gravity  represents  the  effect  of 
material  beneath  the  station,  but  part  of  it  may  be  due  to 
irregularities  in  terrain  (upper  part  of  the  Bouguer  slab)  away 
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from  the  station. 
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d.  Terrain  Effect:  Topographic  relief  around  the  station  has 
a  neqative  effect  on  the  grav i tat ional  force  at  the  station.  A 
nearby  hill  has  upward  gravitational  pull  and  a  nearby  valley 
contributes  less  downward  attraction  than  a  nearby  material 
would  have.  Therefore,  the  corrections  are  always  positive. 
Corrections  are  made  to  the  SBA  when  the  terrain  effects  were 
0.1  milligal  or  larger.  Terrain  corrected  Bouguer  values  are 
called  the  Complete  Bouguer  Anomaly  (CBA).  When  the  CBA  is 
obtained,  the  reduction  of  gravity  at  individual  measurement 
points  (stations)  is  complete. 

A 1.5  INTERPRETATION 

To  interpret  the  gravity  data,  the  portion  of  the  CBA  that  might 
be  caused  by  the  light-weight,  basin-fill  material  must  be  sep¬ 
arated  from  that  caused  by  the  heavier  bedrock  material  which 
forms  the  surrounding  mountains  and  presumably  the  basin  floor. 
The  first  step  is  to  create  a  regional  field.  A  regional  field 
is  an  estimation  of  the  values  the  CBA  would  have  had  if  the 
light-weight  sediments  (the  anomaly)  had  not  been  there.  Since 
the  valley-fill  sediments  are  absent  at  the  stations  read  in  the 
mountains,  one  approach  is  to  use  the  CBA  values  at  bedrock  sta¬ 
tions  as  the  basis  for  constructing  a  second  order  polynomial 
surface  to  represent  a  regional  field  over  the  valley. 

Where  there  are  insufficient  bedrock  stations  to  define  a  satis¬ 
factory  regional  trend,  another  approach  is  to  estimate  the  re¬ 
gional  by  the  process  of  upward  continuation  of  the  CBA  field. 
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In  Potential  Theory,  a  field  normal  to  a  surface,  regardless  of 
its  actual  source,  may  be  considered  as  originating  in  an  areal 
distribution  of  mass  on  that  surface.  if  the  field  strength  is 
known  the  surface  density  of  mass  (grams  per  square  centimeter) 
can  be  calculated.  The  observed  gravity  field  at  the  surface  of 
the  earth  approximately  fulfills  the  requirements  of  this 
theory:  thus  the  observed  (Bouguer  anomaly)  field  can  be  used 
to  compute  a  surficial  distribution  of  mass  which  would  repro¬ 
duce  the  field,  and  most  importantly,  account  for  the  gravity 
field  anywhere  above  the  surface  of  observation.  On  this  basis, 
the  Bouguer  anomaly  field  is  readily  "continued"  to  level  sur¬ 
faces  above  the  ground. 

An  important  property  of  such  "upward  continuation"  is  that  the 
resultant  field  (which  can  be  represented  by  a  contour  map), 
with  increasing  altitudes  of  continuation,  changes  more  with 
respect  to  shallow  sources  than  it  does  with  respect  to  deeper 
sources.  The  anomalous  parts  of  the  field  ascribed  to  shallow 
density  distribution  tend  to  vanish  as  the  continuation  is  car¬ 
ried  upward  whereas  the  field  produced  by  deeper  sources  changes 
only  slightly,  so  that  upward  continuations  produce  "regional"- 
type  fields. 

The  difference  between  the  CBA  and  the  regional  field  is  called 
the  "residual"  field  or  residual  anomaly.  The  residual  field  is 
the  interpreter's  estimation  of  the  gravitational  effect  of  the 
geologic  anomaly.  The  zero  value  of  the  residual  anomaly  is  not 
exactly  at  the  rock  outcrop  line  but  at  some  distance  on  the 
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"rock"  side  ( > i  tiu>  contact.  The  reason  tor  this  1  r;  found  in  the 
explanation  of  the  terrain  effect.  There  is  a  component  of 
gravitational  attraction  from  material  which  is  not  directly 
beneath  a  point. 

If  the  "regional"  is  well  chosen,  the  magnitude  of  the  residual 
anomaly  is  a  function  of  the  thickness  of  the  anomalous  (fill) 
material  and  the  density  contrast.  The  density  contrast  is 
the  difference  in  density  between  the  alluvial  and  bedrock 
material.  If  this  contrast  were  known,  an  accurate  calculation  of 
the  thickness  could  be  made.  In  most  cases,  the  densities  are 
not  well  known  and  they  also  vary  within  the  study  area.  In 
these  cases,  it  is  necessary  to  use  typical  densities  for  ma¬ 
terials  similar  to  those  in  the  study  area. 

If  the  selected  average  density  contrast  is  smaller  than  the 
actual  density  contrast,  the  computed  depth  to  bedrock  will  be 
greater  than  the  actual  depth  and  vice-versa.  The  computed 
depth  is  inversely  proportional  to  the  density  contrast.  A  ten 
percent  error  in  density  contrast  produces  a  ten  percent  error 
in  computed  depth.  An  iterative  computer  program  is  used  to 
calculate  a  subsurface  model  which  will  yield  a  gravitational 
field  to  match  (approximately)  the  residual  gravity  anomaly. 

The  second  vertical  derivative  (SVD)  of  gravitational  field  is 
used  to  aid  the  interpreter  in  evaluating  the  subsurface  mass 
distribution.  Once  the  CBA  field  has  been  projected  onto  a  uni¬ 
form  grid  system,  its  SVD  at  the  grid  nodes  is  readily  computed. 
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In  accordance  with  I, a  Place's  Equation  in  Free  Space,  the  neqa- 
tivo  of  the  second  vertical  derivative  is  equal  to  the  sums 
ot  the  second  derivatives  in  the  x-di rect ion  and  in  the  y- 
direction.  The  second  vertical  derivative  is  an  indication  of 
the  curvature  of  the  Bouguer  anomaly  field.  In  particular  the 
zero-value  of  the  SVD  indicates  the  inflection  in  the  field  as 
it  changes  from  "concave-upward"  (algebraically  negative  SVD ;  to 
"convex-upward"  (algebraically  positive  SVD).  In  a  general  way 
the  zero  SVD  falls  on  the  tightest  contours  of  the  field  and 
where  contours  are  nearly  parallel  its  location  can  be  estab¬ 
lished  by  eye.  However,  where  contours  diverge,  converge,  or 
change  direction  this  is  not  always  so  readily  done.  The  zero 
SVD  contour  line  may  be  an  indicator  of  a  line  of  faulting,  the 
pinchout  of  a  stratum,  truncation  of  a  stratum  at  an  unconform¬ 
ity  or  merely  a  marked  change  in  shape  or  in  density  of  a 
geologic  unit. 
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APPENDIX  A2.0 
PINE  VALLEY,  UTAH 
GRAVITY  DATA 
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